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Spin accumulation and spin transfer torques induced by Spin Hall Effect in bi-layer structures 
comprising ferromagnetic and paramagnetic materials are theoretically investigated. The charge and 
spin diffusion equations taking into account spin-flip and spin Hall effect are formulated and solved 
analytically and numerically for in structures with simplified and complex geometry, respectively. It 
is demonstrated that spin torques could be efficiently produced by means of Spin Hall effect which 
may be further enhanced by modifying structure geometry. 

INTRODUCTION 

In [l|, 0] Dyakonov and Perel predicted the existence of Spin Hall Effect (SHE) in paramagnetic metal without 
applying the external magnetic field under influence of spin-orbit interaction. It was found later a lot of similarities 
between the "anomalous" Hall Effect in ferromagnetic metals (FM) and SHE 0, 0|, assuming that in both cases 
the mechanisms of creation of the transversal electric current is the skew scattering or side jump mechanism due to 
the spin-orbit interaction. The more detailed analyses of SHE was presented in [5|-|7| employing the semi-classical 
Boltzmann equation in Q and Keldish formalism in 0, 0| for the calculation of the transport properties of the 
paramagnetic metal taking into account the spin-orbit interaction. It is considered that SHE may be very effective 
tool for the manipulation with spin current and spin accumulation. The most interesting could be the hybrid structure 
consisting of ferromagnetic (FM) metal/paramagnetic (PM) metal with large SHE, if both mechanisms of creation of 
the spin current are involved: SHE and spin polarization by ferromagnetic metal. In this work we investigate spin 
transfer torques produced by the interplay between these mechanisms in FM/NM bi-layer structures by solving charge 
and spin diffusion equations. 



MODEL 



For the calculation both spin accumulation and spin polarized current in such a structure which may have a complex 
geometry, it is necessary to derive spin diffusion equation taking into account both SHE for the paramagnetic metal 
and processes governing the spin transport in ferromagnetic metal. It is necessary also to develop the convenient code 
for the numerical simulations for the system of the complicated geometry in 2D and 3D cases. Following [8- 11 [ where 
the diffusion equation describing spin transport in ferromagnetic multilayered structures where developed and with 
diffusion equation obtained in [7[ and describing SHE we can write down: 
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where cto is the conductivity, f3 is the spin-asymmetry parameter of conductivity, ctsh is the spin Hall conductivity, 
v - density of states, Um = M/M s , where M is magnetization vector in ferromagnetic, Um = m/\rn\, where to is 
spin accumulation vector, index i is a component of vectors to, j m and Um in spin space, l s { - spin diffusion length 
and lj - exchange spin diffusion length. Inside the ferromagnetic metal one have to omit the terms of SH, and in 
paramagnetic metal (3 = 0. In equation (p} we added the last term, which is quadratic on the V<£> as the value of to is 
proportional to V(f, and in J2]) we omitted terms corresponding to the contribution of the anomalous velocity (see (2) 
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in Q). Here we have to mention the articl e fl2| . where it was proven that the large ctsh/oo, experimentally observed 
for Au doped by Fe and Pt impurities fl~3lll4j and for Cu doped by It may be attributed to the resonant electron 
scattering on impurities if take into account spin-orbit interaction. 

Let us consider the system, consisting of two flat layers, one of the paramagnetic metal with SHE and second 
ferromagnetic layer with current in x direction. If to consider the case L x > L y 3> + , where L x and + 
are the lengths of the system in x and z directions, the solution of ^ in the region of x far from (—L x , L x ) may be 
easily found, and expression for ip, (m^ = m^ 3 ' = 0) are the following: 
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where Li and L 2 - are thicknesses of SH- layer and FM-layer respectively. 




FIG. 1: Schematic of Pt/Py bylayer. Sizes are in "nm". 1 
Magnetisation of Py is in xy plane at n/4 angle to x axe. 



Pt layer, 2 - Py layer, 3 - Cu electrodes. Current is along x axe. 



The spin accumulation produces the effective field H c g acting on the magnetisation of the ferromagnetic layer, 
which value is equal H c h = J s dj /J-b, where J s d is s-d exchange integral. It is important to notice that this field is 
proportional to drop of voltage and not to current density, as in the case of Oersted field created by the current. So 
if one choose as a source of SHE dirty paramagnetic metal only due to its higher resistance the value of the induced 
by SHE effective field for the constant current density will increase, and besides that the following conclusion of 12 1 
the value of ctsh/cto may increase as well. Another interesting conclusion, following from expression for the potential 
(p, is that SHE produce drop of voltage in z direction perpendicular to the current. If there is no ferromagnetic layer 
above SHE structure the drop of voltage is symmetric (cx cosh[z/Z s f]) and quadratic on the applied voltage V, but in 
presence of FM layer this drop has linear on V part and is finite across the thickness of the system. To investigate 
the influence of the edge of layers on spin accumulation in the system of finite size we have solved equation of ([3|) 
numerically using Comsol Multiphysics for several artificial system, consisting of SHE substrate and ferromagnetic 
layers or dots situated on the top of this substrate. 
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FIG. 2: Effective fields induced by Pt SHE in Py near the interface. Adopted values of the parameters: a pt — 0.005(f2 ■ nm) 
I™ = Wnm, cr^ = 0.1cr p ', a Py = 0.0022(0 ■ nm) -1 , l Py = 6nm, /3 = 0.7, lj = Iran, current density j = 10 7 A/cm 2 . 




FIG. 3: Schematic of Pt/Py bylayer. Sizes are in "nm". 1 - Pt layer, 2 - Py layer, 3 - Cu electrodes. Current is along x axe. 
Magnetisation of Py is in xy plane at 7r/4 angle to x axe or along z axe. 



RESULTS 

In FigjT] we give the schematic of Pt/Py bilayer, which was experimentally investigated in [l6j |. In figH] the 
dependence on x coordinate of the induced by SHE fields -ffsHE acting on the magnetization inside of Py layer and 
near the Pt/Py interface is shown. It is clear that SHE spin current produce all tree components of the field, meanwhile 
Oersted field Hj produced by the current in given geometry has only y-componcnt, and Hj = 4 Oe. The z-componcnt 
of the SHE field produces torque in the plane of the Py layer, and the ratio H^tm/Hj lies within the interval 1.3-^3.7, 
where 1.3 is the average value over the volume of Py layer while 3.7 represents this ratio near the interface, what is 
close to the value recalculated using value of S/A = 0.63 in [ll| which gives Hsim/Hj = \J 1 + 47rM c ff / iJ g S/A = 1.9 
in notations of 16]. 

To investigate the possibility to manipulate with magnetization of small ferromagnetic dots under influence of 
SHE current we have calculated the field induced by SHE in the small ferromagnetic dot situated on the surface of 
paramagnetic metallic layer with SHE as shown in Fig. [3l The direction of magnetization of the ferromagnetic dot 
forms the angle 7r/4 with cc-axes, or it is parallel or antiparallel to y axes. In the case M \\ y only j/-component of 
induced fields not zero and its average value does not change under inversion of magnetization direction along y axes. 
This field is much stronger than Oersted field which is about 5 Oe. Besides that field does not produce any torque but 
defines the most energetically favorable direction of magnetization. When this direction is not collinear with y-axes, 
the SHE produces not only y-component of spin accumulation but all three xyz components of spin accumulation due 
to the precession of spin accumulation vector in exchange field of ferromagnet. 

In fig[5]the dependence on x coordinate of three components of spin transfer torque, produced by spin accumulation 
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FIG. 4: Effective fields in Co/Pt structure for the direction of magnetisation parallel and antiparallel to y axe. Only H y 
component survive. a Co = 0.005(f2 • nm)' 1 , = lOnm. 
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FIG. 5: Effective torques for the case of Um = (cos 7r/4; sin 7r/4; 0) near the interface. Averaged values over volume of Co are 
< T x >— 15 Oe, < T x >— — 15 0e, < T z >= 3 0e. 



due to the SHE are shown. The torque acting on magnetisation is defined as M x 7?she > an d in LLG equation 
it has to be multiplied by the gyromagnetic ratio 7. For the case U m = (cos7r/4; sin7r/4; 0) the xy part of torque 
vector T = (T x ; T y ; 0) is lying in xy plane perpendicular to U m and may by considered like an additional damping or 
antidamping term in LLG equation. The torque due to Oersted field has only z-component. 

For the case of M \\ Oz, T x and T y torque components value close to the interface are shown in figlU These 
torques try to reorient magnetization of FM layer onto xy plane and T y component may be considered like damping 
or antidamping term. We checked the results taking the spin Hall conductivity equal to zero and obtained that all 
torques besides one produced by Oersted field of the current vanish. 
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FIG. 6: Effective torques near the interface in case of M \\ Oz. Averaged values over volume of Co are < T x >= — 4 0e, 
< T y >= -30 Oe. 



We have shown that SHE may represent a powerful tool for the manipulation with magnetization of small ferromag- 
netic metal dots situated on the surface of the thin paramagnetic metal layer with large value of SHE conductivity. 
Especially important feature of spin transfer torques created by spin accumulation due to influence of SHE is that 
these torques have component similar to damping or antidamping spin torque produced by current in non-collinear 
magnetic multilayers. 
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